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The origin of the iron lines in NGC 7213 
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Abstract. The analysis of a simultaneous XMM-A'ew;o«/BeppoSAX plus three previous BeppoSAX observations revealed 
that NGC 7213 is a rather peculiar Seyfert 1. No significant Compton reflection component was observed, while an iron line 
complex, best explained in terms of three narrow lines, is clearly apparent in the data. Due to the absence of the reflection 
component, the neutral iron line is likely not produced in a Compton-thick material, but current data do not allow to choose 
between a Compton-thin torus and the BLR. The two ionized iron lines from FeXXV and FeXXVI may be produced in a 
photoionized gas with a column density of a few 10 23 cirT 2 , in analogy with the case of NGC 5506. 
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1. Introduction 

Iron Ka emission lines are common features in Seyfert Is. The 
profile of the line provides fundamental information on their 
origin. If it is produced in the innermost regions of the accre- 
tion disk, kinematic and relativistic effects contribute to forge 
a bro ad, asymmetric, double-peaked profile (see IFabian et alJ 
2000j for a review). A systematic study of the iron hne profiles 
in more than 20 Seyfert Is observed by ASCA suggested that 
most of them ar e likely to have broad lines dNandra et all 19971 
lRevnoldsll997l) . 

Rather surprisingly, this scenario is not being confirmed 
by XMM-Newton and Chandra. Relativistic ally broa dened 
lines a re detected in only a few sour ces ( Nandra et al. [ 1 19991 
IWilms et alJl200ll iTurner et alJl2002t IFabian et alJl2002l> . and 
often the profiles and/or time behaviour are not what expected 
from simple models. The reason for this lack of observed rela- 
tivistic lines may reside in the physical properties, e.g. the ion- 
ization state, of the accretion disk itself. 

On the other hand, a narrow component of the iron line 
is almost invariably present i n Chandra high energy grat ings 



and XMM CCD spectra (e g.lKaspi et alj|2001 



Yaqoob et al 



Gondoin et al 



200 1| iReeves et alJ l200lt IPounds et alJl200ll 
200 it iMatt et alJl200l1) 7The line, typically unresolved at the 
energy resolution of the CCDs, must be produced far from 
the nu cleus, either in th e torus envisaged in the Unification 
model dAntonuccill 1 9931) or in the Broad Line Region. In the 
former case, and if the torus is Compton-thick, a Compton re- 
flection component should also be present, which can be best 
studied when simultaneous XMM/BeppoSAX observations are 
available, in order to exploit the unprecedented sensitivity of 
XMM at the iron line energy and the hard X-ray coverage of 



Table 1. The log of all analysed observations and exposure 
times. 



Date 


Mission 


Instrument 


Texp (ks) 


05/30/1999 


BeppoSAX 


MECS 


47 






PDS 


56 


11/20/1999 


BeppoSAX 


MECS 


15 






PDS 


34 


05/03/2000 


BeppoSAX 


MECS 


46 






PDS 


31 


05/27/2001 


BeppoSAX 


MECS 


61 






PDS 


38 


05/29/2001 


XMM-Newton 


EPIC-PN 


30 



the BeppoSAX PDS instrument. On the other hand, if the line 
is produced in the BLR, a much fainter Compton reflection 
component is expected, and the intrinsic width of the iron line 
should be the same as that of the optical broad lines. 

The nuclear activity of NGC 72 13 (z=0.005977) wa s 
discovered by the HEAO A-2 satellite jMarshall et alJll979l) . 
It was then included in the IPiccinotti et alJ i 19821) sample 
and classified as a Seyfert 1 because of the presence of 
broad optical e mission lines and a strong blue continuum 
dPhillipslll979l iFilippenko & Halpernl \l9S4 . In this Letter, 
we will discuss the X-ray spectrum of the source, analysing a 
simultaneous XMM-Afewfo«/BeppoSAX observation and three 
further BeppoSAX observations, yet unpublished. 
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We will assume Hq = 70 km/s/Mpc throughout the paper. 
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2. Observations and data reduction 

2. 1. XMM-Newton 

NGC 7213 was observed by XMM-Newton on May 29 2001 
(see Table 03 with the imaging CCD ca meras, the EPIC-p n 
JstruderetalJl200ll) and the EPIC-MOS dTurner et alJl200ll) . 
adopting the Medium filter and operating in Prime Small 
Window and Prime Partial W2 respectively. The selected subar- 
ray allows a pileup free spectrum for the pn, since the observed 
count rate is much lower than the maximum defined for a 1% 
pileup (see Table 3 of the XMM-Newton Users' Handbook). 
On the other hand, the MOS count rates are only slightly lower 
than the maximum value: a more careful inspection of the data 
with the S AS tool EPATPLOT, together with a comparison with 
the pn spectra, suggests that the MOS data are probably af- 
fected by pileup. Therefore, we will not use the MOS data 
in this Letter. The data were reduced with SAS 5.4.1. X-ray 
events corresponding to pattern 0-4 were used for the pn and an 
extraction radius of 40" was chosen for spectra and lightcurves. 
Spectra were analysed with XSPEC 11.2.0. 

2.2. BeppoSAX 

BeppoSAX observed the source four times, from 1999 to 
2001 (see Table 03- The last observation was performed 
between May 27 and May 30 2001, thus being simultaneous 
to the XMM-Newton one. Data were downloaded from the 
ASDC archive 1 and analysed with XSPEC 11.2.0. Since the 
softer part of the spectrum is beyond the scope of this work, 
LECS data will not be treated in this Letter. A normalization 
factor of 86 was adopted between the PDS and the MECS 
jFiore et alJ Il999l) . appropriate for PDS spectra extracted 
with fixed rise time threshold. When dealing with the fit of 
the simultaneous XMM-Afewton/BeppoSAX observation, we 
adopted a normalization factor of 0.888 between the EPIC-pn 
and the PDS, which was obtained normalizing the pn and the 
MECS spectra in a combined fit. 

In the following, errors are at the 90% confidence level for 
one interesting parameter (Ay 2 = 2.71), where not otherwise 
stated. 



3. Data analysis 

3. 1. Flux and spectral variability 

The source showed a long-term variation in the 2-10 keV flux 
from 2 to 4x1 0~ 11 erg cirT 2 s (corresponding to a lumi- 
nosity ranging from 1.5 to 3xl0 42 erg s" 1 ) when comparing 
the BeppoSAX observations. Short-term variations always less 
than 5% around the mean value were also observed within sin- 
gle BeppoSAX observations. No significant differences in the 
continuum slope, that could be associated with the flux vari- 
ability, are apparent in the data. 



3.2. Spectral analysis 

We first tried to fit the pn/PDS simultaneous observation in the 
energy band 0.5-220 keV, with a simple power law with an 
exponential cut off, absorbed by the Galactic column density 
2.04 x 10 2l) cm 2 , as found with the Heasarc nH tool 2 . The 
result is an unacceptable value of the x 2 (1500/261 d.o.f.). A 
soft excess and complex absorption features are present at en- 
ergies below ^ 2.5 keV, together with residuals around the iron 
line energy. On the other hand, no significant excess is left at 
higher energies. In this Letter, we will focus on the iron line 
properties and the Compton reflection component: therefore, 
we will ignore the spectrum below 2.5 keV in the following 
analysis, where not otherwise stated. 

The residuals around the iron line clearly suggest that it is 
complex (Fig. 03- If fitted with a single gaussian line, the cen- 
troid energy is 6.39^j jjj keV and the intrinsic width cr is < 90 
eV. The inclusion of a second gaussian line, with E — 
keV and cr < 250 eV, improves the fit at the 99% confidence 
level, according to F-test. A third narrow line (6.94^^) may 
also be present at the 96% confidence level. The fit is now per- 
fectly acceptable, with three narrow (cr fixed to 1 eV) gaussian 
lines (173/178 d.o.f.). 

The presence of a cold iron line, if produced by Compton- 
thick matter, should come along with a significant Compton 
reflection com ponent. Therefor e, we decide d to adopt the 
PEXRAV model (Magdz iarz & Zdziarskill995l) . which includes 
this component in terms of the solid angle R=Q/27r subtended 
by the reflecting material. The fit does not improve at all 
(173/177 d.o.f.), and the Compton reflection component is 
completely undetect ed (see T able Et- We further tried an ion- 
ized reflector model (Ballant vne et al.l200lh . to reproduce both 
the weak Compton hump and the ionized component of the iron 
line: the resulting^- 2 is significantly worse than for the previous 
model (204/179 d.o.f.). 

To verify if the lack of detection of the Compton reflection 
is affected by excluding the softer part of the data, we extended 
the best fit model to the broad band spectrum (0.5-220 keV). As 
already mentioned, residuals in excess to the model are clearly 
present below 2.5 keV. After allowing for a new minimization 
of the^ 2 , the fit is still unacceptable (y 2 =539/256 d.o.f.), and 
absorption features are left in the softer part of the spectrum. 
The addition of a warm absorber (model ABSORI in XSPEC) 
improves marginally the fit and only an upper limit on the col- 
umn density of the gas is found. On the other hand, if the resid- 
uals below 2.5 keV are treated like a genuine soft excess, in- 
cluding a black body to the original model, a good^- 2 (287/252 
d.o.f.) is achieved. In this case, the other spectral parameters, 
notably F and R, remain unaltered with respect to the fit re- 
ported in Table |2 

To look for the presence of the reflection component in the 
past history of the source, we decided to analyse all the avail- 
able BeppoSAX observations (see Table Q3. We adopted the 
best fit model used for the simultaneous observation, again ig- 
noring data below 2.5 keV. There is no significant detection of 
Compton reflection component in any of the observations, the 



1 http://www.asdc.asi.it/ 



2 http://heasarc.gsfc.nasa.gov/cgi-bin/Tools/w3nh/w3nh.pl 
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6.2 6.7 
Energy (keV) 



Fig. 1. EPIC pn residuals around the iron line when no gaussian 
lines are added to the model (see text for details). The energies 
are not corrected for the redshift of the source. 




Energy (keV) 



Fig. 2. Best fit model and residuals for the pn/PDS spectrum 
(see text for details). 



Table 2. Best fit parameters for the pn/PDS simultaneous ob- 
servation. Th absorbing column density was fixed to that of our 
Galaxy, the intrinsic width of the lines was kept at the value 
of 1 eV because they are unresolved and a value of i = 30° 
was adopted for the unconstrained inclination angle (see text 
for details). 



(keV) 



r 

e c 

R 

Ej (keV) 

EWi (eV) 

Fi (10 5 ph cm 2 i 

E 2 (keV) 
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1.69 ±0.04 
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6.39 ± 0.01 
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2 5 +io 
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0.6±0.3 
6 94+ 05 

: 14 
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173/177 
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0.5: 




value of R being always consistent with zero. Since no variation 
in the continuum slope is observed between the observations, 
we performed a combined fit of all of them, with all the param- 
eters linked but the normalization factors needed to account for 
the different flux levels. In Fig.[5]we show the F-R contour plots 
for this combined fit: the 90% (for two parameters of interest) 
upper limit on R is ^ 0.2. This result is unique among bright 
Seyfert Is observed by Bep poSAX (e.g.lPerola et al.ll2002l) . A 
similar case was claimed by Weaver (2001): a very tight upper 
limit on the Compton reflection component in MCG-2-58-22 
was found, based on ASCA and RXTE data, together with a sig- 
nificant iron line. However, results from a BeppoSAX/XMM- 
Newton si multaneous observation suggest a higher value of R 
feianchi etall2003hl) . 

Finally, we checked for a possible evidence of variability on 
the 6.4 keV line, leaving free the line normalizations between 
the different BeppoSAX observations in the combined fit: no 
significant improvement in the x 1 is found. 

4. Discussion 

The 6.4 keV line is clearly too narrow to be produced in the 
innermost region of the accretion disc, which would require 
a significant broadening of the feature (see Fig. |4j. Its ori- 



1.64 1.66 1.68 1.7 1.72 1.74 1.76 1.78 

r 

Fig. 3. T-R contour plot for the combined BeppoSAX observa- 
tions. The curves refer to Ay 2 =2.30, 4.61, 9.21, corresponding 
to confidence levels of 68, 90 and 99% for two interesting pa- 
rameters. 

gin must be placed in a material farther away from the BH, 
either Compton-thick, such as the 'torus', or Compton-thin, 
like the Broad Line Region. However, the notable absence of 
a Compton reflection component, as reported in the previous 
section, argues against an origin from Compton-thick material. 
In fact, a line with EW^ 80 eV, as observed by XM M-Newton 
would require a reflection com ponent with R =* 0.5 dMatt et alJ 
199 it iGeorge & Fabian 111 99 ll) . The only way to reconcile the 
data with the presence of a Compton-thick material would be 
assuming a torus with a small covering fa ctor and an iron over- 
abundance of about 3-5 (Ma tt et alll9971) . at odds with typica l 
values of about 1 found in Seyfert galaxies jPerolaetalJ2 002). 

Th erefore, the line s hould originate in a Compton-thin ma- 
terial. iMatt et alJ (120031) have shown that an iron line with 
EW^ 80 eV can be produced by a face-on Compton-thin torus 
with Nh — 2 x 10 23 crrT 2 . In this case, the associated reflec- 
tion component would be much less prominent than for the 
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6.34 6.36 6.38 6.4 6.42 6.44 



E (keV) 

Fig. 4. cr-E contour plot for the 6.4 keV iron line in the simul- 
taneous BeppoSAX/XMM-Afewfcw observation, when the line 
width is left free to vary. The curves refer to A^ 2 =2.30, 4.61, 
9.21, corresponding to confidence levels of 68, 90 and 99% for 
two interesting parameters. 



Compton-thick material. A similar conclusion can be drawn 
if the line is produced in the BLR. To discriminate between 
the two possibilities, it is in principle possible to rely on the 
intrinsic width of the feature, which is significantly higher if 
the line is produced in the BLR, where the gas keplerian ve- 
locity is larger. Unfortunately, the 90% (for two parameters 
of interest) EPIC-pn upper limit to the line FWHM is 8 000 
km s _1 (see Fig.|4}, which is similar to the broader co mponent 
of the , indeed rather complex, Ha line ( Filipp enko & Hal pern 
11984 . thus providing no further information on the origin of 
the line. Only a Chandra gratings observation could resolve 
the line width or put a tighter upper limit. No further informa- 
tions can be derived from line variability. As reported in the 
previous section, the line flux is consistent with being constant 
among all the analysed observations, but the error bars are very 
large. 

The observed centroid energies of the other two narrow 
lines indicate that they are produced by Fe XXV and Fe XXVI. 
The two ionized lines must therefore originate in a medium 
different from that o f the neutral line, in analogy with th e 
case of NGC 5506 llMatt et alJl200lt IBianchi et alJl2003ah . 
Bianchi & Matt] J2002I) have shown that a gas photoionized by 
a power law continuum can produce both lines, with EWs of a 
few tens of eV with respect to the total continuum: assuming 
a covering factor f=l and a solar iron abundance, the required 
column density is around 10 23 cirT 2 . A photoionized gas with 
such a column density would produce a continuum component 
with the same spectral shape of the primary emission and a flux 
of around 15%. However, differently from NGC 5506, the pri- 
mary continuum in NGC 7213 is not absorbed, preventing to 
disentangle the two components. A study of the variability of 
the ionized iron lines in response to the variation of the nuclear 
continuum would be useful in this case, but the BeppoSAX ob- 
servations can provide only upper limits to the fluxes of these 
lines, making the investigation impossible. 



5. Conclusions 

NGC 7213 appears unique among Seyfert Is observed by 
BeppoSAX, because it lacks a significant Compton reflection 
component in all of the four available observations performed 
by the Italian-Dutch satellite. The upper limit on R for the com- 
bined fit of the BeppoSAX observations is 0.2. Moreover, the 
EPIC pn spectrum, simultaneous with the last BeppoSAX ob- 
servation, reveals the presence of an iron line complex, best ex- 
plained in terms of a narrow neutral iron line with EW^80 eV 
and two ionized lines from Fe XXV and Fe XXVI, with EW of 
a few tens of eV. Because of the absence of a Compton reflec- 
tion component, the neutral iron line is likely produced in the 
BLR or a Compton-thin torus, even if current data do not allow 
to make a choice between these two alternatives. The ionized 
lines may be produced in a photoionized gas with a column 
density of a few 10 23 cm -2 , but other solutions could not be 
excluded. 
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